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Abstract. - We have calculated the electron density maps of the ripple phase of dimyris- 
toylphosphatidylcholine (DMPC) and palmitoyl-oleoyl phosphatidylcholine (POPC) multibi- 
layers at different temperatures and fixed relative humidity. Our analysis establishes, for the 
first time, the existence of an average tilt of the hydrocarbon chains of the lipid molecules along 
the direction of the ripple wave vector, which we believe is responsible for the occurrence of 
asymmetric ripples in these systems. 



Lipids self-assemble in water to form a variety of lamellar phases 0, 0|. The ripple or 
Ppi phase characterized by a one-dimensional height modulation of the bilayers is seen in 
some phospholipids under high hydration In the phase diagrams of these systems, 

it is sandwiched between the high-temperature L a phase and the low-temperature L^/ phase 
[|[ H . In the L a phase, the hydrocarbon chains of the lipid molecules are in a molten state 
and the in-plane ordering is liquid-like. On the other hand, in the Lpi phase, the chains are 
predominantly in the all- trans conformation and are tilted with respect to the layer normal. 
The chains are also ordered in the plane of the bilayer, but the exact nature of the degree of 
this ordering is yet to be determined ^ || . 

Experimental studies on the ripple phase have established many of its structural features. 
Almost all x-ray studies show asymmetric ripples corresponding to an oblique unit cell of 
the rippled bilayers Jl], ^, [|, ||; though there have been some reports of symmetric ripples 
corresponding to a rectangular unit cell . The latter have been shown to be metastable 

structures in some systems Q. It is well established that only lipids that have a Lp> phase 
at lower temperatures exhibit the ripple phase, indicating the importance of the tilt of the 
chains in the formation of the ripples Q . Determination of the chain tilt in the ripple phase 
is of utmost relevance as it is the key structural feature hitherto unknown. But it has not 
been possible to obtain it directly from x-ray diffraction patterns as the chain reflections are 
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rather diffuse, probabiy due to the presence of disordered chains, ie., chains that are not 
entirely in the trans conformation. NMR and diffusion experiments also detect a population 
of disordered chains in this phase ||, [l(| . Hence detailed information about chain tilt has to 
be deduced from the electron density map of the system, calculated from x-ray diffraction data. 

The lack of knowledge of its structural features has hindered the formulation of a satisfactory 
theory of the ripple phase; none of the current theories is consistent with all the experimental 
observations. The most striking disagreement concerns the occurrence of asymmetric ripples. 
It has been proposed that molecular chirality is responsible for such ripples [ |TT| , but experi- 
ments indicate otherwise [Q, |l^| . 

Recently, Sun et al. jl3| have calculated the electron density map of the ripple phase of 
DMPC using the x-ray data of Wack and Webb g. They find that the ripples have a 
saw-tooth shape, with the bilayer thickness in the minor arm being much smaller than that 
in the major arm. Further, the electron density in the headgroup region along the major arm 
is much higher than that along the minor arm. This led them to hypothesize that the chain 
organization in the major arm is like in the Lpt phase and that in the minor arm is like in 
the L a phase. However, this conjecture is not supported by other experiments. For example, 
self-diffusion in the ripple phase is found to be highly anisotropic, with a fast component that 
is 4-5 orders of magnitude faster than the slow component Jl(| ; but the fast component itself 
is about 2-3 orders of magnitude smaller than that in the L a phase. Thus the authors of Ref. 
10 conclude that although the intramolecular hydrocarbon chain disorder may be substantial 
in the fast bands, the intermolecular order in this region is not like that in the L a phase. 

In view of this discrepancy, we have calculated the electron density maps of the ripple phase 
of DMPC and POPC. In addition to the x-ray data of Ref. 5, we have used data from oriented 
samples at different temperatures and fixed relative humidity. We find that the ripples in both 
these systems have a saw-tooth shape, with the ratio of the lengths of the two arms essentially 
independent of temperature. If the molecules in the short arm were in the L a phase, the 
length of this arm would be expected to increase as the L a phase is approached from below. 
This is contrary to what we see. Further, the difference in the bilayer thickness and the 
electron density in the two arms can be largely accounted for in terms of an average tilt of 
the chains along the direction of rippling, which we believe is responsible for the occurrence 
of asymmetric ripples in these systems. These results are clearly important for a satisfactory 
theoretical description of this phase. 

We have adopted the modeling and least squares fitting procedure developed by Sun et 
al. jn| to calculate the electron density maps. The unit cell parameters of the two-dimensional 
oblique lattice are the two vectors a and b, and the angle 7. In terms of the ripple wave length 
A and the lamellar spacing d, the two lattice vectors can be expressed as: a = d cot 7 x + 
d z, and b = A x. Here x is the direction of the ripple wave vector and z is the direction of 
the average layer normal (see Fig.l). A, d, and 7 are directly measured from the diffraction 
pattern. The electron density within the unit cell, /?(x,z), is described as the convolution of a 
ripple contour function C(x,z) and the transbilayer electron density profile T^(x,z). C(x,z) = 
<5(z-u(x)), where u(x) describes the ripple profile and is taken to have the form of a saw-tooth 
with peak-to-peak amplitude A. Ai is the projection of the longer arm of the saw-tooth on the 
x-axis. T^(x,z) gives the electron density at any point (x,z) along a straight line, which makes 
an angle ip with the z-axis. The electron density in the methylene region of the bilayer is close 
to that of water and is taken as zero. 
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We have used three models for T^(x,z), two of which are equivalent to the SDF and M1G 
models of Ref. 13. In model I, it is taken as consisting of two delta functions with positive 
coefficients pu, corresponding to the headgroup regions separated by a distance L, and a 
central delta function with negative coefficient of magnitude pm, corresponding to the methyl 
region. The six adjustable parameters in model I are: A, Ax, ip, Ph/ ' Pm, L and a normalizing 
factor. In model II, the delta functions representing the head and methyl groups are replaced 
with Gaussians of width 07, and <r m , respectively. The electron density in the minor arm is 
allowed to be different by a factor fx from that in the major arm. The region where the two 
arms meet is modeled as a wall with an electron density differing by a factor from the 
rest of the arm. The wall thickness in this model is fixed at a small value. Thus there are 



10 adjustable parameters in this model. It is using these two models that Sun et al. |13 
find that the bilayer thickness along the local layer normal is different in the two arms of the 
ripple. But this result is built into these models as the parameter L, which is the thickness of 
the bilayer along a direction that makes an angle ip with the z-axis, is taken to be the same 
in the two arms. Therefore, in model III, we remove this constraint and allow L as well as 
ah, a rn and ph/pm to be different in the two arms of the ripple. Further, the wall between 
the two arms is taken to have a variable width w. This model has 15 adjustable parameters. 
Minimization is done by iterative least squares fitting with respect to six variables at a time 
& 

The structure factors at the observed (h,k) values are calculated using each of the above 
models. These are then compared with the observed structure factors and a chi-square value is 
obtained, which is subsequently minimized by varying the adjustable parameters in the model. 
The phase of each of the Bragg reflections is obtained from the structure factors calculated from 
the converged model. These calculated phases are combined with the observed magnitudes 
of the structure factors and inverse Fourier transformed to get the electron density map of 
the system. We have used the x-ray diffraction data of Wack and Webb [0 from powder 
samples of DMPC as well as our data from oriented films of Z-DMPC, dZ-DMPC and POPC. 
Details of the experimental procedure are discussed elsewhere A . Relevant geometric intensity 
corrections were applied to the data from oriented samples, but absorption corrections could 
not be applied as the sample thicknesses are not accurately known. We have confirmed by 
assuming reasonable values of the sample thickness, that these corrections do not significantly 
change the calculated electron density profiles. However, in the absence of these corrections, 
we are unable to analyze these data using models II and III, due to the lack of convergence 
resulting from the large number of adjustable parameters in these models. 

The converged values of the different parameters in the three models are given in table 1, 
along with the crystallographic R factor. The precision of these parameters is about 0.1 %. 
The data from Ref. 5 were used in these analyses. In accordance with the results of Ref. 
13, we find that all the phases, except those of the three relatively weak (0,k) reflections, 
obtained using model II are the same as those obtained from model I. These phases are also 
the same as those obtained from the SDF and M1G models in Ref. 13. Model III gives only 
a marginally better fit to the data, the chi-square being only about 2% lower. The phases of 
all the reflections other than the (0,k) reflections are the same as those found from models I 
and II. For the (0,k) reflections this model gives a combination of phases that is different from 
the ones obtained from models I and II. The values of L in the two arms are almost the same 
and equal to that obtained from model II. This is also true of the other parameters which 
are allowed to be different in the two arms. However, model III gives a slightly higher value 
of 0.7 for fx- As discussed below, this factor can be accounted for in terms of the chain tilt, 
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without resorting to the assumption of a L Q -like organization in the minor arm. The low y 2 
and R values for models II and III, and the absence of any physically unacceptable features 
in the electron density map (see Fig. 2) indicate that these models closely represent the true 
structure of the system. 

The electron density map of the ripple phase of DMPC, calculated with the data of Ref. 
5 is shown in Fig. 2. The ripples clearly have a saw-tooth shape, with an offset between the 
two leaves of the bilayer. The simplest explanation for this offset is an average tilt of the 
chains along the rippling direction; such an offset cannot be expected if the tilt were in a plane 
normal to the rippling direction. The tilt angle ip is found to be approximately equal to ("y-w). 
Further confirmation of the existence of an average tilt along this direction comes from the 
fact that the value of L is almost equal in the two arms and is comparable to twice the length 
of a fully stretched DMPC molecule. If it is assumed that the chains are tilted at an angle 
tp with the z-axis, their tilt with respect to the local layer normal can be calculated from the 
shape of the ripple. Using the values of the structural parameters given in table I, the tilt 
angle with respect to the local layer normal turns out to be 1.6° and 34.5° in the longer and 
shorter arms, respectively. The tilt in the short arm is comparable to that found in the Lpi 
phase. Since the area per molecule is inversely proportional to the cosine of this angle, a value 
of 0.82 is obtained for f\. This is in very good agreement with the value of 0.77 obtained from 
the map for the ratio of the average electron densities in the headgroup region of the longer 
and shorter arms. Thus an average tilt of the chains along the rippling direction provides a 
consistent explanation for many features of the electron density map. This means that to a 
good approximation the height modulation of the bilayers along the x-axis can be described 
as arising from a relative sliding movement of neighboring chains, with all the chains lying 
the x-z plane and tilted by a constant angle ip with respect to the z-axis. The existence of an 
average chain tilt along the rippling direction breaks the reflection symmetry of the bilayer in 
the plane normal to it and hence can be expected to be responsible for the asymmetric ripples 
seen in this system. 

We have also calculated the electron density maps of /-DMPC, dZ-DMPC and POPC at 
different temperatures in the ripple phase, using data from oriented films. The structural 
features of the ripples are found to be similar to those obtained from the data of Ref. 5. 
The maps of the chiral and racemic DMPC samples were identical, indicating the lack of 
influence of molecular chirality on the ripple structure p2| . The temperature dependence 
of the structural parameters of the ripples in DMPC are found to be very weak, as in the 
case of dipalmitoyl phosphatidylcholine (DPPC) |l5). Contrary to what is observed in freeze 
fracture experiments fig] , we find that the ripple shape has Fourier components higher than 
the second. Further, we do not find a significant temperature dependence of the amplitude of 
the ripples in contrast to what is reported in Ref. 17. 

The electron density map of the ripple phase of POPC is shown in Fig. 3. The ripple 
shape is very similar to that of DMPC. It also has a saw-tooth shape and an offset between 
the monolayers indicating an average tilt in the direction of rippling. In POPC, the angle 7 is 
much larger than in DMPC, whereas the wave length and layer spacing are comparable. Unlike 
those of DMPC, the structural features of the ripple phase of POPC vary significantly with 
temperature, as shown in table 2. In the absence of absorption corrections, the fits are not as 
good as in the case of DMPC data of Ref. 5. As mentioned earlier, we find that the electron 
density profiles are insensitive to these corrections. Hence the values of the last two parameters 
quoted in the table are the ones estimated from the electron density maps. The layer spacing 
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decreases slowly and 7 increases steadily as temperature is increased. The ripple wavelength 
first decreases and then suddenly increases to a large value just below the transition. These 
trends are very similar to those seen in DPPC GEL but in POPC the temperature dependence 
is very much pronounced. The amplitude, except near the L Q transition, is about half that 
of DMPC. In both DMPC and POPC the ratio of the projected lengths of the major and 
minor arms is about 2 and is essentially insensitive to temperature. This observation further 
supports the view that the chain organization in the minor arm is not like that in the L a phase. 

All the freeze fracture studies of the ripple phase show ripples oriented over micrometer-sized 
regions [jl], |l6|, [l?], [l8| . Since the chain tilt is locked to the rippling direction, this implies 
long range order of the tilt direction. These experiments also show ripples oriented only 
along three directions, each at an angle of approximately 120° from the other two, indicating 
a six-fold symmetry in the underlying bilayer structure. Thus the in-plane ordering of the 
molecules in the bilayer is at least hexatic. 

In conclusion, we have calculated the electron density maps of the ripple phase of DMPC 
and POPC. The shape of the ripples in these two systems are very similar, with both of them 
exhibiting asymmetric ripples. We have been able to establish the existence of an average 
chain tilt in the direction of rippling, which is probably responsible for the asymmetric ripples 
seen in these systems. 

We thank Y. Hatwalne, K. Usha, and J. F. Nagle for discussions and HKL Research Inc. 
for the use of their software. 
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Table I. - The values of the structural parameters obtained from the three models. 



Parameter 


Model I 


Model II 


Model III 


A(A) 


19.7 


18.7 


19.1 


A X (A) 


106.2 


101.0 


101.5 


j, 


5.7° 


9.9° 


9.1° 


L(A) 


41.0 


37.0 


37.2,37.0* 


Phi P-m 


1.1 


0.9 


1.1,1.1* 


fl 




0.6 


0.7 


h 




0.9 


0.65 


Oh 




4.9 


4.6,4.7* 


C m 




13.3 


9.4,10.0* 


W(A) 




2.0(fbced) 


9.2 


x 2 


692.5 


217.6 


213.6 


R 


0.172 


0.083 


0.087 



* indicates the values in the shorter arm 



Table II. - Temperature variation of the structural parameters of the ripple phase of POPC at 
75%RH. 



T(°C) 


d(A) 


7° 


A(A) 


A (A) 


Ai 
(A-Ai) 


13.0 


58.3±0.1 


116±1 


200±2 


10±1 


2.0±0.1 


13.5 


58.0 


119 


170 


9 


1.8 


14.5 


57.3 


124 


143 


5 


2.6 


15.0 


56.4 


133 


266 


5 


1.9 



V X 




Schematic diagram showing the structural 
parameters of the ripple phase 



Fig. 1. - Schematic showing the structural parameters of the ripple phase. 
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x(A) 

Fig. 2. - The electron density map of the ripple phase of DMPC obtained using the data from Ref.13. 
T=18.2°C and the volume fraction of water is 0.263. A=141.7A and 7=98.4°. The positive (negative) 
contours are represented by solid (dotted) lines. The regions with positive electron density correspond 
to the head groups. Note that the thickness of the bilayers is about 40 A, whereas that of the water 
region is about 20 A. 




Fig. 3. - The electron density map of the ripple phase of POPC obtained using x-ray data from 
oriented samples. T=13.5°C and RH=75%. 



